JLAB-THY-11-8 



Deeply Virtual Exclusive Processes and Generalized 
Parton Distributions 

Charles E. Hyde^'^, Michel Guidal^, and Anatoly V. Radyushkin^'^'^ 

^ Laboratoire de Physique Corpusculaire, Universite Blaise Pascal, 63177 Aubiere, FRANCE 

^ Department of Physics, Old Dominion University, Norfolk VA, 23529, USA 

^ Institut de Physique Nucleaire ORSAY, 91406 Orsay, France 

* Thomas Jefferson National Accelerator Facility, Newport News, VA, 23606, USA 

^ Bogoliubov Laboratory of Theoretical Physics, 141980 Dubna, Russia 

Abstract. 

The goal of the comprehensive program in Deeply Virtual Exclusive Scattering at Jefferson 
Laboratory is to create transverse spatial images of quarks and gluons as a function of their 
longitudinal momentum fraction in the proton, the neutron, and in nuclei. These functions are 
the Generalized Parton Distributions (GPDs) of the target nucleus. Cross section measurements 
of the Deeply Virtual Compton Scattering (DVCS) reaction ep — > ep7 in Hall A support 
the QCD factorization of the scattering amplitude for > 2 GeV^. Quasi-free neutron- 
DVCS measurements on the Deuteron indicate sensitivity to the quark angular momentum sum 
rule. Fully exclusive H(e, e'pj) measurements have been made in a wide kinematic range in 
CLAS with polarized beam, and with both unpolarized and longitudinally polarized targets. 
Existing models are qualitatively consistent with the JLab data, but there is a clear need for 
less constrained models. Deeply virtual vector meson production is studied in CLAS. The 12 
GeV upgrade will be essential for for these channels. The p and ui channels reactions offer the 
prospect of flavor sensitivity to the quark GPDs, while the (/)-production channel is dominated 
by the gluon distribution. 



1. Deeply Virtual Exclusive Scattering 

In the past decade, Deep Exclusive Scattering (DES) has emerged as a powerful new probe of 
the partonic structure of the nucleon, hadrons and nuclei. These are reactions of the type: 

e + Nucleon — )• e + Nucleon + 7 

— )• e + Nucleon + meson (1) 

The eN — )• eN^ reaction is the coherent sum of the Bethe-Heitler and virtual Compton 
amplitudes, as illustrated in Fig. [TJ The electron scattering kinematics of deeply virtual 
processes corresponds to the Deep Inelastic Scattering (DIS), or Bjorken limit of inclusive 
electron scattering, with large and W'^ above the resonance region. In addition, the 
forward exclusive limit is defined kinematically by —t < Aq^j-, < 1 GeV^. Thus the deeply virtual 
Compton scattering (DVCS) amplitude, 7*A'^ — >• 7A'^ is an "off-forward" generalization of the 
forward Compton amplitude which defines the DIS cross section via the optical theorem. 

The intense interest in DVCS started after the article by Ji [T], linking DVCS to the total 
contribution of quarks to the proton spin. It was found [H O [3l H] that, in analogy with DIS, 



in the limit of large and small t, amplitudes of DVCS and deeply virtual meson production 
(DVMP) can be expressed in a power series of with the power determined by twist 

(dimension minus spin) of each operator in the expansion. Detailed proofs of factorization for 
DVCS and DVMP were given in O [6l [71 [8]. This factorization is depicted in Fig. [2] for the 
DVCS and DVMP amplitudes. In these processes, the leading power term of the amplitude is 
the convolution of the perturbative kernel with a new class of non-local bi-linear twist-2 quark 
(or gluon) operators, called Generalized Parton Distributions (GPDs) [Il[2]. These were first 
described by Miiller et al. [9]. In the case of deeply virtual meson production, the hard kernel 
is also convoluted with the meson Distribution Amplitude (DA). 

1.1. Generalized Parton Distributions 

The Generalized Parton Distributions (GPDs) parameterize Fourier transforms of nucleon 
matrix elements of bilinear quark (and gluon) operators separated by a light-like interval 
= f9]. The kinematics are commonly defined in terms of symmetric variables (see Figs. 
Dandig) : 




Figure 1. Lowest order QED amplitude for the ep — )• ep^ reaction. The momentum four- 
vectors of all external particles are labeled at left. The net four- momentum transfer to the 
proton is = {q — q')fi = {p' — p)fi- In the virtual Compton scattering (VCS) amplitude, the 
(spacelike) virtuality of the incident photon is = —q^ = —{k — k')'^. In the Bethe-Heitler 
(BH) amplitude, the virtuality of the incident photon is — = —t. Standard (e, e') invariants 
are Sg = (A; + p)'^, x^j = Q'^/i^q ■ p) and W'^ = {q +p)'^. 




Figure 2. Factorization of the 7*j> — >• 7p DVCS amplitude and the 7*p — )■ MA'" deep virtual 
meson production amplitude in the Bjorken limit of large and —t <C Q^. All permutations 
of the photon and gluon vertices should be included in the amplitudes. The labels on the quark 
lines are the light-cone momentum fractions, relative to P"^ = {p + p')^/2. 



The generalized Bjorken variable has the same form with respect to the symmetrized variables 
P and q as does xbj with respect to the DIS variables p and q. 

It is convenient to use a reference frame in which P'^ has only time- and z-components, both 
positive. We define light-cone vectors 

n^^ = [1, 0, 0, -1] /{V2P+) , r = [1, 0, 0, 1] P+/{V2) (3) 

Then in the forward limit of either DVCS or deeply virtual production of a light meson, —2,^ is 
the "-|-" fraction of both the momentum transfer to the target and the virtual photon: 

A+ = A • n w -2^P+ g . n = g+ (4) 

The quark GPDs H and E are the nucleon helicity conserving and helicity-flip matrix elements 
of the vector operator containing 7 • n = 7"*". Suppressing the QCD scale dependence and the 
Wilson-line gauge link one can write the flavor-/ dependent GPDs as (see, e.g., |10j): 



/ 
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Hfix, t) 7 • « + Ef{x, e, t) ^n„(T"^A^ 



U{p,s) , (5) 



where the U{p, s) are the nucleon spinors. The factorization proofs demonstrate that the initial 
and final -|- momenta of the active parton are (x zizS,)P~^. The GPDs H, E are defined similarly 
as the matrix elements of the axial operator containing n ■ 775 = 7~'"75: 



dz P ^i^p+, )j-n-f5^f{z )\p,s) 



2-K 

= U{p\s') 



~ ~ n • A 

Hf{x, t)n ■ 775 + Ef{x, ^, *)^^75 



Uip,s). (6) 



With the convention that positive and negative momentum fractions refer to quarks and anti- 
quarks, respectively, we observe the following kinematic regions of the GPDs (Fig. [2]): x > ^ > 0: 
the initial and final partons are quarks; x < — ^ < 0: the initial and final partons are anti-quarks; 
\x\ < a qq pair is exchanged in the t-channel. This identification is refiected in the QCD 
evolution equations of the GPDs. For |x| > ^, the evolution of GPDs is similar to the DGLAP 
evolution of the forward parton distributions, whereas for \x\ < ^, the GPDs evolve according 
to the ERBL equations of a meson DA [2]. 

The GPDs combine the momentum fraction information of the forward parton distributions 
(PDFs) of DIS with the transverse spatial information of the elastic electro-weak form factors. 
The forward limits of the helicity conserving GPDs are 



^/(x,0,0)] (q^^x)9{x)-qf{-x)9{-x) 
Hf{x, 0, 0) J \Aqf{x)9{x) + Ag^(-x)0(-x) 



(7) 



where g/(x) and qf{x) are the fiavor-/ dependent quark and anti-quark momentum-fraction 
distributions; and Aqj(x) and Agj(x) are the quark and anti-quark helicity distributions. There 
are no specific analogous constraints on the forward limits of E and E. 

The first moments of the GPDs are equal to the corresponding elastic form factors. 




where Fij{—t), F2j{—t) are flavor-/ components of tlie Dirac and Pauli form factors of the 
proton, defined with positive arguments in the space-like regime. Similarly, 



\ J [gpji-t), 

where qaj and gpj are the flavor components of the axial and pseudoscalar form factors of the 
proton. 

The .^-independent integrals of ^ and ^ are examples of the polynomiality condition 
required by Lorentz invariance. Specifically, the moment of a GPD is polynomial in even 
powers of ^, with maximal power < {N + 1) for H, E and maximal power < N for iJ, E. The 
second moment of the GPD sum H + E leads to the important angular momentum sum rule of 
JiP: 



lim / dxx[Hf{x,i,t)+Ef{x,i,t)] =2Jf, 



(10) 



where 2Jf is the fraction of the spin of the proton carried by quarks of flavor /, including 
both spin and orbital angular momentum. More generally, at non-zero t, the individual second 
moments are form factors (i.e., Fourier transforms of spatial distributions) of the nucleon's 
energy- momentum tensor [l]. 

fi 4 

/ dxxHf{x,C,t) = M2.f{t) + -fdijit) 

Jo -3 

/■I 4 

dxxEf{x,^,t) = [2Jj{t)-M2j{t)]--edij{t) . (11) 

The forward limit M2j(0) is the ordinary momentum sum rule 

M2j{0)= [ xdx[qf{x) + qf{x)] . (12) 
Jo 

The term dij{t) is Fourier conjugate to the spatial distribution of the time-averaged shear stress 
on quarks in the nucleon [111 [T2] . 

The GPDs themselves, and not just the moments, provide unique spatial information about 
partons in the hadronic target. The parton impact parameter b is Fourier conjugate to A_|_. 
The Fourier transform of Hf{x,0, A'^) determines a positive-definite probability distribution of 
quarks of flavor / as a function of longitudinal momentum fraction x and spatial coordinate b in 
the transverse plane \13\ 114]. As a consequence, the Dirac form factors Ff'^ are the 2D Fourier 
transforms of the charge densities of the proton and neutron in impact parameter space. A 
recent analysis reveals the presence of a negative charge density at the heart of the neutron [15]. 
In the context of the Double Distribution models of the GPDs (see section II. 3p . this can be 
understood by the excess of down quarks over up quarks in the neutron at large x. Similarly 
to H, the combination of Hf and at ^ = determines the spatial density of quarks in 
a transversely polarized proton [16]. These distributions strongly break azimuthal symmetry 
about the longitudinal axis. In particular, the centroid of the up and down distributions are 
displaced in opposite directions, as required by the fact that the Lorentz boost of a magnetic 
dipole produces an electric dipole field. For 7^ 0, the Fourier transform of the GPDs determines 
overlap matrix elements for partons of initial and final impact parameter b/(l it^) with respect 
to the center-of-momentum of the initial and final proton [T^. In the particular x = case, the 



variable r, Fourier conjugate to A_|_, is the transverse separation of the active parton from the 
center-of-momentum of the spectator partons |18| . 

The experimental program to establish the domain of factorization in DES and to extract 
GPDs promises new insight into the quark-gluon structure of hadrons. The GPDs offer for the 
first time a probe of the rich correlations between spatial and momentum degrees of freedom of 
quarks and gluons in hadronic systems. 



1.2. Scattering Amplitude and Observables 

The ep — )• epj cross section has the form (Fig. [T]) 
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(13) 



where ^%ycs ~ ^x^-^M"^ / Q"^ and X is the BH-VCS interference term. The pure Bethe-Heitler 
term |7bhP is exactly calculable in terms of the nucleon form factors [19 1 120 ^ [21]. The full VCS 
amplitude is expressed as 



rvcs(e±) =n(fc',A)7;.n(fc,A)M/fM-4 . 



(14) 



The general VCS hadronic tensor H has 12 independent terms. In the leading order twist-2 
approximation, H reduces to just four terms [T] 



-'-'LO,twist2 



i i-gn± U{p') [{n . j)nit t) + ^n,a-^ AxSiC, t)] U{p) 
- {en^ U{p') [{n ■ J^5)m, t) + (75n • A)£{C, t)] U{p) , 

where the transverse tensors are defined as 



(15) 



(16) 



The Compton form factors (CFF) % . . . in (jl5p are defined by the integration over the quark 
loop in Fig. [2) 
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(17) 



Thus the imaginary parts of the CFFs are proportional to the GPDs at the point x = it^. 
Complete expressions for the VCS hadronic tensor to twist 3 accuracy are given in e.q.[2H [T2]. 

The importance of the azimuthal distributions of the DVCS and interference terms, both for 
testing factorization and extracting constraints on the GPDs was first pointed out by Diehl et 
aZ.[22]. The azimuthal distribution of the DVCS and interference terms has the general form |20j 
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(18) 
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The azimuthal angle ^-y-y of the hadronic q' (8) p' plane relative to the electron scattering plane 



is defined such that sin( 



> when (k A k') • q' > and 



■'77 



when the final photon is on 



the beam side of q. 

The leading order twist-2 DVCS amplitude couples only to transverse photons. Consequently, 
there is a specific twist-hierarchy to the terms in (|18|19p |20j : 

(i) cjj'^'^s, 4 are twist-2; 

(ii) (c, s)^ are twist-2 (transverse VCS in interference with longitudinal BH amplitudes); 

(iii) (c, 5)5^^^^ are twist-3 (LT electroproduction interference terms); 

(iv) {c,s)^^'~'^, are bilinear combinations of ordinary twist-2 and gluon transversity terms; 

(v) (c, 3)2 are twist-3; 

(vi) (c, ,3)3 are linear in the twist-2 gluon transversity terms. 

For example, the c, terms for unpolarized and longitudinally polarized targets are respectively: 
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(20) 



(21) 



where A is the target polarization and tq = — A^/(4M^) [20]. The EM form factors Fi, F2, Gm 
are evaluated at —t and the Compton form factors %,£... are evaluated at (^, t). In particular, 
we expect that on the proton the F\H and FiH terms will dominate the unpolarized and 
longitudinal target polarization observables of (j20p and (j2ip . respectively. On the other hand, 
on the neutron both Fi and gi{x) = H(x,0,0) are small. We anticipate a greater sensitivity 
to E from the interference term on an unpolarized neutron and to the combination of Ti and 
iS on a longitudinally polarized neutron. Transverse target observables depend on different 
combinations of CFFs [20]. Dynamic and kinematic twist-3 terms in the scattering amplitude 
cause kinematically suppressed twist-3 terms to mix into the "twist-2" observables listed above 
[23^ I24j. In addition, any "twist-2" observable will naturally contain contributions of all higher 
even- twist in a power series in 1/Q^. The precise twist content of DES observables must be 
determined via a Q^-dependent analysis at fixed {^,t). A complete analysis must also include 
the logarithmic Q^-dependence from QCD evolution. 



1.3. Models 

Several models of the GPDs exist. To varying degrees, they incorporate the theoretical and 
empirical constraints on the GPDs. In the valence region, the most widely used models are based 
on the Double Distribution (DD) ansatz proposed by Radyushkin [25]. Detailed versions of 
this model are presented by Vanderhaeghen, Guichon, Guidal (VGG) [26] and Goeke, Polyakov, 
Vanderhaeghen [21]. The Double Distributions re-parameterize the (x,^) dependence of the 
GPDs in terms of the momentum fractions /? and a, of and A"*", respectively. Thus the 
initial and final parton + components of momentum are /3P'^ =F (1 i a)A"'"/2. The H, E, and 



H Double Distributions are parameterized as: 



GPDf^DD{x,^,t) = dl3 da6{x - (3 - aOFf{P,a,t) 

J-i J-i+m 

t^fqfim-PTf/Af, (22) 
I A5^(/3); 



F/(/3,a,0) = /i(/3,a) < 



for H, E, and H, respectively. In (|22|) qf and Aqj are the ordinary and helicity-dependent 
parton distribution functions of flavor / and is the flavor anomalous magnetic moment of the 
proton. The normalization of E is such that Aj = J — P)^^ qf{P)- The profile function h is 
commonly parameterized as 

HB a) - r(26 + 2) 



For 6=1, this form reduces at /3 = to an asymptotic meson DA <I>(a) = 3(1 — a^)/4, with 
support — 1 < a < 1. This connection is suggested by the ERBL evolution equations for |x| < ^. 
In general, the exponent 6 is a free parameter and for 6—7-00, the GPD is ^-independent. 

The DD form of (j22p ensures that the polynomiality conditions are automatically satisfied. 
However, it was pointed out by Polyakov and Weiss that for H and an additional "L'-term" 
must be included [27] to produce the highest ^^+^ power for moment. This term, which only 
has support in the ERBL region, is an iso-singlet and enters with opposite sign to H and E: 

Hfix, t) = Hf,nD{x, ^,t) + e{C- \x\) j^D {x/t t) , 

Ef{x,tt) = Ef^DD{x,^,t)-9i^-\x\)^D{x/tt) . (24) 

iV f 

In practice, the D-term has been taken as an expansion in odd Gcgenbauer polynomials, with 
the first few terms fitted to a Chiral Soliton model calculation [28. ,2T] . The t dependence is 
introduced into the model via a Regge inspired ansatz |21j : 

Ff{P,a,t) =/i(/3,a)gj(/3)|/3|-"'* . (25) 

As found in [29], the high-t form factor is dominated by the contribution in the DD at large /?, 
for which the simple Regge forms must be modified to describe the data. A fit to all of the 
nucleon form factors data was obtained with the ansatz 

Ff{f3,a,t) =/i(/3,a)gy(/3)|/3|-(i-/^)"'*. (26) 

Aside from the choice of 6-parameter for each flavor and GPD, variants of the model exist with 
different choices of including the valence or valence plus sea contributions to E and H. The GPD 
E is generally parameterized separately, as the pion-pole in the t-channel. In this framework, E 
is directly related to the pion form factor [2T| I30j. 

The family of models sketched above, and generically labeled "VGG" is qualitatively 
successful in describing the DVCS data. However, the model is highly constrained and does 
not have the full degrees of freedom of the GPDs. More general parameterizations will be 
needed as the data improves in precision and covers both a broader kinematic range and a more 
complete set of spin and flavor observables. 



Another approach [311 [32] is to construct valence generahzed parton distributions as an 
overlap of light-cone wave functions. However, a model involving the lowest Fock state 
component only produces GPDs vanishing at the border points x = and in the whole 
central region < x < ^ \31\ I33j . It was shown [H^ that inclusion of the higher Fock state 
components gives GPDs that are nonzero in the central region and at the border points. In 
particular, one may assume that overlap of the lowest Fock state components gives model GPDs 
at a low normalization point Qq ~ 300 MeV, and then evolve them to hard scales Q ^ 1 GeV: the 
evolution will induce nonzero values for GPDs in the central region. Originally [35], the evolution 
approach was used to build a model for the gluon GPD, assuming that for a low normalization 
point Qq it coincides with the usual ("forward") gluon density, H3{x,^;Qq) = xG{x,Qo). In 
Ref. |36j this ansatz was also applied for quark distributions in an attempt to describe HERA 
DVCS data at low xbj, for which predictions based on the double distribution ansatz are too 
large in magnitude. More recently, the "Dual Parameterization" (DP) framework developed by 
Polyakov and collaborators \37\ [38] was used to address this issue. In this approach, GPDs are 
expanded in terms of the partial waves exchanged in the t-channel. It was expected that for 
low XBj of the HERA DVCS data, the expansion may be truncated to the first "forward-like" 
functions [38]. However, a detailed analysis [39| demonstrated that the minimal model of the 
dual parameterization significantly (by a factor of 4) overestimates the HERA data. The relation 
between the dual parameterization approach and the double distribution ansatz was investigated 
in Ref. |40] . where it was shown that GPDs built from DD-based models with 6 = 1 in Eq. (j23p 
and small may be reproduced just by the first term of the dual parameterization expansion, i.e., 
the minimal DP- and DD-based models give similar results for DVCS at small xbj, and both give 
a rather large value ~ 1.8 for the ratio -R^(C) = -ff^(C) 0/^(0 of singlet quark distributions for 
small ^, while experimental data favor the value close to 1. In a model developed by D. Mueller 
and collaborators [H] it is possible to keep the value of i?^(0 "fiexible", i.e. to adjust it to 
describe the data. The "flexibility" may be achieved also in the dual parametrization approach, 
if one adds the second "forward-like" function. 

2. Initial DVCS Experiments 

The richness of physical information in the GPDs has sparked an intense experimental effort. 
The HI, HERMES, and CLAS collaborations published the first evidence for the DVCS reaction 
in 2001. 

2.1. DVCS at HERA 

The HI [42l US] E] and ZEUS [45l [46] collaborations measured the p{e, e'j)X cross section, 
integrated over (p^^. The exclusive p{e,e''y)p channel is enhanced over p{e,e''y)N* channels by 
vetoing on forward detectors [44l|46]. In ZEUS, a subset oi p{e,e'jp) events were tagged in a 
forward tracker [46]. The HERA data cover a wide kinematic range at low xbj, with central 
values of and W from 8 to 85 GeV^ and 45 to 130 GeV, respectively. 

The HERMES collaboration measured the (/>^^-distribution of the relative beam- 
helicity asymmetry in the H(e, e'7)X reaction at average kinematics (Q^, xbj, t) = 
(2.6 GeV^, 0.11, -0.27 GeV^) [47]. The FWHM of the Mj. distribution was ^ 1 GeV^ therefore 
covering the majority of the resonance region. However, at low xbj and —t, model estimates 
indicate that the exclusive H(e, e'7)p channel is dominant [47] . The HERMES collaboration has 
recently measured the beam-charge asymmetry |48] . transversely polarized target asymmetries 
|49] . longitudinally polarized target asymmetries [50], and a more extensive set of beam spin 
asymmetries |51] . The final 2006-2007 HERMES run utilized a new recoil detector |52] . to 
establish exclusivity via H(e, e'7p) triple coincidence [53] . 



2.2. Initial CLAS DVCS Data 

The JLab CLAS Collaboration first measured the relative beam-helicity asymmetry in the 
H(e, e'p)x reaction with 4.25 GeV incident electrons [5l]. The exclusive photon was detected 
in only a fraction of the acceptance, due to the limited small angle acceptance of the standard 
CLAS calorimeter. The distribution is shown in Fig. [3] (left). The position and width of 
the exclusive H(e, e'p)'y event distribution was constrained to fit a subsample of H(e, e'pj) data 
at small (e'7) opening angle, such that the events are dominated by the BH process. Similarly, 
exclusive vr*^ event distributions were constrained to a subset of H(e, e'p^yy) events from vr^ decay. 
Thus the H(e, e'p)x events in the exclusive region were fitted with two gaussians, for the x = 7 
and X = TT*^ channels. The widths and positions of these two gaussians are a priori constrained. 
In this way the exclusive H(e, e'p)^ channel was isolated. The resulting DVCS beam-helicity 
asymmetry is shown in Fig. [3] (right). The shaded band is a one-sigma fit of the form 

BSA = asm{(j)jj) + Psm{2cl)^^) (27) 

The a coefficient contains the twist-2 physics. The /3 coefficient contains the twist-3 physics, as 
well as contributions from cos{(j)^^) terms in the unpolarized cross section in the denominator of 
the beam spin asymmetry. The dashed and dotted curves in Fig. [3] are leading twist calculation 
of the VGG model, in the .^-independent (at fixed x) and ^-dependent versions, respectively [26] • 
The solid curve includes an estimate of twist-3 effects [55l[56]. The models, though constrained 
by fundamental principles are still very preliminary. It is remarkable that the data and models 
are in as good agreement as indicated by Fig. [3l Within the VGG model, the largest contribution 
to the beam helicity asymmetry on the proton comes from the -ff(ib^, ^, t) GPD. 

A second CLAS experiment, still with the standard CLAS configuration [57], measured the 
longitudinal target spin asymmetry in the p{e, e'p^) reaction on a polarized NII3 target [58] • In 
order to isolate the exclusive ep — )■ ep7 events from the nuclear continuum, the statistics were 
limited to the triple coincidence (e, e'p'j) events, with the photons detected in the standard CLAS 
calorimeter. The resulting exclusivity spectrum in Fig.|3]shows a 10:1 signal to background ratio. 
The longitudinal target spin asymmetry, averaged over the acceptance, is displayed in Fig. [5l for 
= 1.82 GeV^ (0 = 0.16, and (t) = -0.31 GeV^ The solid curve is a fit of the same form 
as (j27p . The resulting sin(0^^) moments are plotted in Fig. [6l The error bars in Figs. [5] and [6] 
are statistical, with the systematic errors displayed as a band at the bottom. The dashed and 
dotted curves in Figs. [5] and [6] indicate the sensitivity of the longitudinal target spin asymmetry 
to H. 

The initial success of the CLAS DVCS analysis, and the limited small angle acceptance of the 
CLAS calorimeter led to the construction of a small angle "Inner Calorimeter". The ongoing 
dedicated DVCS program in CLAS will be described in section HI 

3. The Hall A DVCS Program at 6 GeV 

The Hah A DVCS program started with experiments E00-110[59] and E03-106[60]. These 
experiments measured, respectively, the cross sections of the II(e, e'^)p and D(e, e'^)pn reactions 
at XBj = 0.36 with an incident beam of 5.75 GeV. In both experiments the scattered electron 
was detected in the standard High Resolution Spectrometer (HRS) [61], and the photon was 
detected in a new 132 element PbF2 calorimeter, subtending ~ 0.1 sr. PbF2 is a pure Cerenkov 
medium, thereby minimizing the hadronic background and delivering the fastest timing pulses. 
All PbF2 channels were readout by a custom 1 GHz digitizer [62 j. based on the ANTARES ARSO 
chip|63j. The luminosity of 1-4 •lO'^^Hz/cm^ per nucleon was unprecedented for open detectors 
in a non-magnetic environment. The halo-free CW beam of CEBAF was essential to this success. 




Figure 3. }l{e,e'p)x analysis from CLAS at 4.25 GeV incident electron energy |54| . 
Left: Missing mass squared distribution. The two Gaussian distributions representing the 
H(e, e'p)j and H(e, e'p)7r^ events are described in the text. The smooth polynomial background 
represents the contribution of processes such as ep — )■ eN*j. Right: DVCS Beam-helicity 
asymmetry. The kinematics are integrated over 

Q2 g [1., 1.75] GeV^ and -t G [0.1, 0.3] GeV^. 
The shaded region is the fit described in the text, in section [2^21 The curves, described in the 
text, are evaluated at the fixed values Q"^ = 1.25 GeV^, xbj = 0.19, and t = — O.lDGeV^. 



3.1. Proton DVCS 

Hall A experiment EOO-110 measured DVCS on the proton at = 1.5, 1.9, and 2.3 GeV^. 
The isolation of the exclusive H(e, e'j)p signal is illustrated in Fig. [71 The helicity dependent 
cross sections as a function of (f)jj in four bins in are displayed in Figs. [8] and [H The latter 
figure also displays the helicity independent cross sections for = 2.3 GeV^. The helicity 
dependent cross sections demonstrate the dominance of the effective twist-2 term of (fT9]) . 
The helicity independent cross sections (Fig. [9]) show significant contributions from the sum of 
the interference and DVCS terms, in addition to the pure BH cross section. Thus the analysis 
of relative asymmetries of the form Aa/a requires the inclusion of the full DVCS terms in 
both the numerator and denominator. The effective "twist-2" interference term QmC^ of (j20p 
is presented in Fig. [TOl The VGG model calculation, described in section 11.31 agrees in slope 
with the data, but lies roughly 30% above the data. Within statistics, the results in Fig. [10] are 
close to Q^-independent in all bins in A^. This provides support to the conjecture that DVCS 
factorization results in leading twist dominance at the same scale of > 2GeV^ as in DIS. 

3.2. Neutron DVCS 

JLab Hall A experiment E03-106 measured the helicity dependent DVCS cross section on 
deuterium, D{e,e'j)X at = 1.9 GeV^ and xbj = 0.36. Within the impulse approximation, 
the cross section is described as the incoherent sum of coherent deuteron and quasi free proton 
and neutron channels: 



D(e, e'7) = d{e, e'^)d + n(e, e'^)n + p(e, e^)p + 



(28) 




Figure 4. Missing mass squared distributions of the (e, e'p)x reaction for (e, e'p7) events 
on a longitudinally polarized NH3 target. Left: Raw distribution; Right: Distribution after 
requiring detected photon within 2° of the predicted direction of an exclusive photon from the 
(e, e'p)7 kinematics. In both Figs, the stars are the corresponding distributions from a C target, 
normalized to the negative region. 




Figure 5. CLAS Longitudinal target spin asymmetry, Aul [SH], for (Q^) = 1.82 GeV^ 
(0 = 0.16, and (t) = -0.31 GeV^. The solid curve is a fit of the form of The dashed 

and dotted curves are from the ^-dependent VGG model with E = E = 0. The dotted curve 
includes only H. The dashed curve includes both H and H. 
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Figure 6. The sm{(j)^^) moments of Ajjl 



Left: Three bins in —t, integrated over ^; 



Right: Three bins in ^, integrated over t. The curves are described in Fig. [5l 
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Figure 7. Missing mass squared distribution of the H(e, e'^)X reaction in JLab Hah A 
experiment EOO-110 [64] ■ The [green] stars are the raw data after accidental subtraction. The 
continuous [black] histogram is the data after subtracting the statistical sample of H(e; e''y)'yX' 
events inferred from the measured H(e, e''iT^)X' sample. The open [red] cross histogram is a 
normalized sample of H(e, e'jp) events. The [magenta] dots are the exclusive simulation. The 
[blue] triangles and squares are obtained by subtracting the last two histograms from the solid 
[black] histogram. 
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Figure 8. Helicity dependent DVCS cross sections A'^a/[dQ'^dxBjdA'^d(p^^] from JLab Hall A 
EOO-110 [M] at = 1.5, and 1.9 GeV^. The mean values of -A^ are, from right to left 0.17, 
0.23, 0.28, and 0.33 GeV^. Each distribution is fitted with the form of (jl9p . with the (effective) 
twist-2 term in red and the complete fit in blue. 
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Figure 9. Helicity dependent (top) and independent (bottom) DVCS cross sections from JLab 
Hall A at = 2.3 GeV^ [6l]. The bins are the same as Fig. [8l as are the curves in the top 



plot. In the bottom plot, the green curve (mostly concave up) is the pure |BH| cross section. 
The magenta curve (mostly concave down) is a fit including the twist-2 and twist-3 terms of 

mm. 
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Figure 10. Imaginary part of the effective interference term C^^^ extracted from tlie helicity 
dependent data of Figs. [8]and[9l 

Meson production cliannels contribute as background. The proton-DVCS contribution is 
calculated by smearing the H(e, e''y)X data by the nucleon momentum distribution in the 
deuteron. This statistical estimate of the proton contribution is subtracted from the data. The 
coherent deuteron and quasi free neutron channels were separated, within statistics, by fitting the 
missing mass distribution with a Monte Carlo simulation of these two channels. This separation 
exploits the fact that for calculated relative to a nucleon target, the quasi free neutron 
spectrum peaks at M'^ w whereas the coherent deuteron peak lies at M'^ + t/2. This analysis 
produced constraints on the neutron and deuteron DVCS^BH interference terms QmfC^j^p] |65j . 
Mazouz et al, [65] fitted the neutron interference signal by varying the parameters of the E 
GPD within the VGG model of [21]. This results in a model dependent constraint on the Ji sum 
rule values of {Jd,Ju), illustrated in Fig. [TTl A similar constraint obtained by the HERMES 
collaboration in DVCS on a transversely polarized proton target is also illustrated in the figure. 
Both of these experimental determinations are essentially constraints on the model at one value 
of XBj, and then the model is integrated over x at fixed ^ to obtain the sum rule estimate. 
Measurements over a more extensive range in ^ with a more complete set of spin observables, 
and models with more degrees of freedom are necessary in order to more fully constrain the sum 
rule with realistic error bars. Lattice QCD calculations, and other phenomenological estimates 
are also illustrated in Fig. [TTJ 

3.3. Future Hall A Program at 6 GeV 

The unpolarized cross sections in Fig. [9] are not fully dominated by the pure BH process. 
The harmonic (p^^ structure of the cross section does not allow the full separation of the 
Ke[DVCS^BH] and |DVCS|2 contributions. These terms can be separated either by the beam 
charge dependence {e.g. [48]) or by measuring the incident energy dependence of the cross 
sections. At fixed Q^, xbj, the DVCS, Interference, and BH terms in the cross section scale 
roughly as : Se : I Experiment E07-007 [TT] will measure the DVCS hehcity 

independent cross sections in the three kinematics of Figs. [8] and [9] at two separate beam energies 
in each kinematics. This will measure the dependence of the separated leading twist-2 and 
twist-3 observables of the 3f?e[DVCStBH] and jDVCSp terms. 

Experiment E08-025 [72] will measure the DVCS cross sections on the deuteron at the same 
(<5^,XBj) value as in E03-106, but at two incident beam energies. Together with an expanded 
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Figure 11. Experimental constraints on the total up and down quark contributions to the 
proton spin. JLab Hall A neutron [65] and HERMES transversely polarized proton [39]. The 
theory/model values are from AHLT [66], QCDSF quenched [67] and unquenched [68] LHPC 
[69] . and Thomas [70] . 

calorimeter to improve the neutral pion subtraction, the two beam energies will allow a more 
complete separation of the DVCS^ and real and imaginary parts of the DVCS-BH interference 
on a quasi- free neutron. This will be an important step towards a full flavor separation of DVCS. 
Both experiments E07-007 and E08-025 are running in Autumn 2010. 

4. The CLAS DVCS Program at 6 GeV 
4.1- Unpolarized Proton Targets 

A new calorimeter of 424 tapered PbW04 crystals was constructed to provide complete 27r 
photon coverage for polar angles from 4.5° to 15°, relative to the beam line. A ~ 5 Tesla 
superconducting solenoid was added at the target, to confine Moeller electrons. The new 
calorimeter is located 60 cm from the target where the solenoid fringe field is still a few 
Tesla. Therefore, the individual crystals were read-out by Avalanche Photo-Diodes. Having 
strongly benefited from the CERN CMS pioneering research and development effort on this 
recent technology, the present CLAS experiment is the first one to use such photodetectors in a 
physics production mode. 

All particles of the ep — t- ep'j reaction final state were detected in CLAS. To ensure exclusivity, 
several cuts were made, a couple of them being illustrated in FigJ121 In spite of these very 
constraining cuts, some contamination from the ep — )• epn^ reaction remained. Indeed, when 
one of the two 7's originating from the decay vr'^ ^ 77 escapes detection and/or has little energy 
(below the 150 MeV threshold of the calorimeter), an event ep — t- 6^7(7) may pass all DVCS cuts 
and become a perfect candidate to be selected as an ep — )• ep^ event. Such "I-7" vr'' background 
can be estimated from Monte-Carlo combined with the actual number of detected "2-7" vr'^'s, 
resulting, depending on the kinematics, in contaminations ranging from 1 to 25%, being 5% in 
average. 

The extensive CLAS dataset of DVCS beam spin asymmetries (BSA) is illustrated in Fig. [T3j 
The blue solid curves are the result of the twist-2 handbag GPD calculation (VGG) including 
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Figure 12. Example of exclusivity cuts (given by the location of the arrow) for the CLAS DVCS 
experiment |73j. Distribution in cone angle O-^y for the ep — ?• epY reaction (left) and in missing 
energy Ex for the ep — t- ep^X reaction (right) before (black-dotted curve) and after (red solid) 
kinematic cuts (including others not displayed here). The thin solid black line represents the 
background from the ep — )• epi:^ events. These distributions are integrated over all kinematics 
variables. 

just the H GPD (26l |29]. The blue dashed curves include the associated twist-3 calculation. 
Although the general trends of the data are reproduced, the model tends to overestimate the 
BSAs. These too large BSAs by the VGG model can come from either an overestimation of 
T-L[^,t) (the dominant factor in the numerator of the BSA) or an underestimation of the CFFs 
associated to the real part of the DVCS amplitude, which contribute predominantly to the 
denominator of the BSA |74| . The dashed black curve (third curve in some panels of FigllSp is 
the result of a Regge model [75] for the DVCS process. As increases, the Regge model drops 
significantly below both the data and the VGG calculations. This experiment was continued in 
2008-2009, which will significantly improve the statistical precision, relative to FigJT3] [76] . 

4-2. Polarized Targets 

In 2009, a new DVCS experiment completed data taking with the longitudinally polarized NH3 
target [77j. Relative to the previous experiment ( |58j and Figs. HHS]), this new experiment 
will improve both the statistics and acceptance by the addition of the new electromagnetic 
calorimeter mentioiied in the previous section. We recall that the target spin asymmetry is 
mostly sensitive to T-L{^,t) and that, therefore, strong constraints on this CFF should arise from 
this experiment, as discussed e.g. in [78] . 

After decades of development, the HD-ice target ran successfully at the BNL-LEGS facility 
in 2005 and 2006. This target has now been transferred to JLab and is being prepared for a 
photo-production run in 2011 |79j . Initial studies of local depolarization by microwaves suggest 
that the spin relaxation times of this target are sufficiently long for the target to operate with 
electron beams in CLAS. An electron beam test is projected for the end of the 2011 photo- 
production run. If successful, a full suite of transverse polarization observables for the DVCS 
process will be feasible in CLAS in 2011. When combined with the cross section and longitudinal 
polarization data, along with the resulting double polarization observables, a full separation of 
the real and imaginary parts of all four Compton Form Factors 'H{^,t), £{^,t), T-l{^,t), and 
£{6,,t) is in principle possible [20] [73]. 
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Figure 13. Beam spin asymmetry sm{(j)^^) moments from the CLAS DVCS experiment |73j . 
Curves are described in the text. 

4-3. Nuclear Targets 

GPDs are also defined for nuclei |8m I81| . One can study effects similar to the EMC effect 
observed for standard inclusive parton distributions functions (PDF) where the PDF of a nucleus 
is not simply the sum of the individual nucleon PDFs. A pioneering experiment [H2\ [55] of 
coherent DVCS on a He target ran with the CLAS detector in 2010. ^He IS a very good 
starting case study as it is dense enough to generate nuclear medium effects, many microscopic 
calculations for its nuclear structure and dynamics exist and, as a global spin-0 object, at 
leading- twist, there is only one GPD. This ^He-DVCS experiment detected the scattered electron 
in CLAS, the final state photon with the PbW04 and standard calorimeters mentioned in 
the previous sections and the recoil nucleus with a radial time-projection chamber [84]. The 
(;^-distribution of the coherent DVCS BSA, up to twist-3 corrections, can yield the real and 
imaginary parts of the Compton Form Factor of the coherent GPD. 

5. Deeply Virtual Meson Production 

GPDs are in principle also accessible through exclusive meson electroproduction (see Fig. [2]). 
With respect to the DVCS process, a few features are proper to the meson channels: 

• The factorization holds only for the longitudinal part of the amplitude which implies to 
separate, experimentally, the transverse and longitudinal parts of the cross section. For 



pseudo-scalar mesons, this can be done through a Rosenbluth separation. For vector mesons, 
this separation can be carried out, relying on the s-channel helicity conservation (SCHC) 
concept, by measuring the angular distribution of the vector meson decay products. 

• In comparison to the DVCS handbag diagram, there is now a perturbative gluon exchange. 
This suggests that factorization will be obtained at a higher scale in exclusive meson 
production than in DVCS. 

• Besides the GPDs, there is another non-perturbative object entering the meson handbag 
diagram, the meson distribution amplitude (DA). It is usually taken as the asymptotic DA 
but it potentially adds a further unknown in the process. 

• As a positive point, the meson channels have the advantage of filtering certain GPDs: the 
vector meson channels are sensitive, at leading twist, only to the H and E GPDs while 
the pseudoscalar mesons are sensitive only to the H and E GPDs. Deeply virtual meson 
production also offers a flavor filter of the GPDs. For example, and uj electroproduction 
are sensitive to different combinations of the up- and down-quark GPDs. 

Exclusive vr'' electroproduction results have been published from CLAS [85] and Hall A |86] . The 
Hall C results on exclusive tt"'" production are discussed in the "Transition to Perturbative QCD" 
chapter of this volume [87]. In this section, we focus on exclusive vector meson production. 

5.1. The channel 

Deeply virtual electroproduction of the p^ was studied by the CLAS collaboration at incident 
energies of 4.2 GeV [88] and 5.75 GeV [89]. To select the e~p — >• e~pp^ vr+vr" channel, 
the scattered electron, the recoil proton and the tt^ were detected. A cut on the missing 



resonant e~p — )■ e~pTT~^TT~ (physical) background, arising for instance from processes such as 
e~p — )• e~TT~A^~^ pTT~^ . These "background" channels led to uncertainties of the order of 
20% to 25% on the extracted cross sections. 

The separation of the longitudinal and transverse parts of the cross section was carried out, as 
mentioned earlier, by studying the angular distribution of the decay pions in the center of mass of 
the TT~^TT~ system. At the same time, by the analysis of various azimuthal angular distributions, 
SCHC was verified experimentally at the ~ 20% level. The longitudinal part of the 7*p — )• pp^ 
cross section, which can in principle lend itself to a GPD interpretation, is displayed in Fig. [T4l 
along with the world data. 

The cross sections clearly exhibit two different behaviors as a function of At low W, the 
cross sections decrease as W increases (xbj decreases) and then begin to rise slowly for > 10 
GeV. These two kinematic regimes can be identified, simply speaking, with regimes of t-channel 
exchange of Reggeon or qq exchange in the former case and of Pomeron or 2-gluon exchange 
in the latter case. The results of the calculations of the JML [93] model, based on Reggeon 
exchange and hadronic degrees of freedom, and of the VGG [26j and GK \90\ [9T] models based 
on GPDs and on the handbag diagram of Fig. [2] are shown in FigJ14[ At lower W values, where 
the new CLAS data lie, it is striking that both the GK and VGG models fail to reproduce the 
data even though they are very successful at large W, even at ~ 2.3 GeV^. In the high- 
W (low XBj) region, the gluon GPD calculations already contain large higher- twist effects in 
the form of intrinsic k± effects. The question then arises whether the higher twist effects have 
a different nature in the region dominated by quark GPDs (low W), or whether the double 
distribution based GPD models are missing an essential contribution. Ideas for such "missing" 
contribution in the D-teim of the GPDs are speculated in [89^ 192] . leading to the thick solid 
curve in Fig. [TH 





Figure 14. World data (for W < 50 GeV) for the reduced cross sections j^^p — )• pp^ as a 
function of W for constant bins (/ib). The dashed curve (GK) [90l [91] and the thin sohd 
curve (VGG) [26] are GPD calculations. The thick solid curve is the VGG calculation with 
the addition of the D-term inspired contribution [89l [92] . The dot-dashed curve is the Regge 
JML [93] calculation. The 5.75 GeV CLAS, 4.2 GeV CLAS, CORNELL, HERMES, FermiLab 
and ZEUS data are respectively from refs. [89], [88], [9l], [95], [96], and [97j . 



5.2. The CO channel 

The w channel was studied in CLAS by detecting the e~p — )■ e~pTT^X and e~p — )■ e^pTr+Tr"^ 
topologies [98]. The former is advantageous to determine total cross sections with high statistics 
and the latter is necessary to measure the distribution of the decay products of the uj to 
separate the longitudinal and transverse parts of the cross section if SCHC is verified. However, 
one important result of this experiment was that many SCHC-violating spin density matrix 
elements were measured to be significantly different from in a; electroproduction. Therefore, 
the longitudinal and transverse parts of the cross sections were not separated. Also, the angular 
analysis revealed the importance of unnatural parity exchange in the t-channel, such as vr" 
exchange. 

In terms of quantum numbers, the vr*^ exchange contribution can be identified to the E GPD. 
In the framework of the JML model jl04j . f-channel vr'' exchange is a major contributor to the 
cross section. The suggested importance of vr^'-exchange is in apparent contradiction with the 
theoretical prediction that, at sufficiently large Q^, exclusive vector meson production should 
be mostly longitudinal and sensitive only to H and E. Therefore, in order to study the GPD 
formalism in u production, it is essential to experimentally isolate the purely longitudinal cross 
section, via a Rosenbluth separation. The VGG calculation of cjl, shown in Fig ll5l lies well 
below the unseparated data. This suggests that a precision extraction of ai via a Rosenbluth 
separation will be a difficult experimental challenge. 
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Figure 15. Left: Total (unseparated) cross section for the reaction 7*p — t- pa;, for 
{W)=2.l GeV (top) and {W)=2.8 GeV (bottom), as functions of Q^. The dotted curve is 
the JML model for the total cross section ax + sctl and the solid (dashed) curves shows the 
result of the JML (VGG) calculation for eaL- Data are from CLAS (full circles) [98], DESY 
(empty diamonds) [99] and Cornell (empty circles) [94] . Right: Longitudinal 'y*p — ?• jx^ cross 
section as a function of W at Q'^=3.8 GeV^. Data from CLAS (open circle) [TOU], HERMES 
(solid circle) |101| . ZEUS (open triangles) |102j and HI (solid square) [103] . (Courtesy of P. 
KroU and S. Goloskokov). 

5.3. The (f) channel 

The e~p — )• e~p(f) K^K~ reaction was identified in CLAS by detecting the scattered electron, 
the recoil proton and the positive kaon and cutting around the missing mass of a kaon [lOOj . 
Relying on the SCHC concept, which was experimentally verified to hold in this channel, the 
longitudinal/transverse separation of the cross section was carried out. Fig [TSl shows the resulting 
longitudinal total cross section at = 3.8 GeV^, along with higher energy HERMES and HERA 
data at comparable Q^. 

Exclusive (p electroproduction on the proton can be interpreted in terms of the handbag 
diagram with gluonic GPDs. Fig [15] shows the result of such calculation in the framework of 
the GK model [90^ I91j. The very good agreement between this GPD calculation and the data 
gives confidence in the way higher twists corrections are handled, i.e. by taking into account 
the intrinsic transverse momentum dependence of the partons in the handbag calculation. 

This set of three experiments delivered the largest ever dataset on vector meson production in 
the large valence region. Although conclusions for the meson channels are more challenging 
than for DVCS, there may be the possibility to interpret the and cp channels in terms of the 
handbag diagram, though with large higher-twist corrections and possibly modifications of the 
Double-Distribution based GPD parametrisations. The higher data from JLab at 12 GeV, 
as well as a global analysis including the larger DVCS dataset anticipated in the coming years 
should greatly clarify the role of factorization in deep virtual vector meson production. 



6. Outlook 

6.1. Jefferson Lab at 12 GeV 

The JLab 12 GeV project offers an unprecedented frontier of intensity and precision for the 
study of deep exchisive scattering. The design luminosity of the upgraded CLAS12 detector 
IS 1035/8 cm , with a large phase space acceptance for simultaneous detection of DVCS and 
deeply virtual meson production channels. At this luminosity, the Hall B dynamic nuclear 
polarization NH3 target will achieve a longitudinal proton polarization of 80%. The Hall A 
and Hall C spectrometers will allow dedicated studies at luminosities > lO^'^/scm^ for neutral 
channels 7, vr" at low t and up to 4 • lO'^'^/scm^ for charged channels vr^, K^. Specific 12 
GeV experiments on hydrogen are approved in Hall A for DVCS (E12-06-114), in Hall B for 
DVCS (E12-06-119), and deep virtual 7r°, 77 production (E12-06-108), and in Hall C for deep 
virtual 7r+ production (E12-06-101, E12-07-105). Detailed descriptions of these experiments 
are available on the Hall A, B, and C web pages at www.jlab.org. The projected kinematic 
range of the DVCS programs in Hall A and B is illustrated in FigllGl With CLAS12, additional 
studies are in progress for measurements of deep virtual vector meson production, neutron DVCS 
via 'D{e,e'jn)p (LOI-09-001), coherent deuteron DVCS (PR-06-015) and DVCS on transversely 
polarized protons. 
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Figure 16. Projected kinematic bins at JLab 12 GeV; Left: CLAS12 kinematics for DVCS 
and DVMP on unpolarized H2 and longitudinally polarized NH3 targets. The colors and density 
are proportional to the relative count rates. Right: Hall A kinematics for DVCS and vr'^ 
electroproduction. Beam time is adjusted for roughly equal counts in all bins. 

6.2. Beyond 12 GeV 

The COMPASS experiment at CERN proposes to measure DVCS in high energy muon scattering 
at low via triple coincidence H(/2^,^^,7p) detection |105j . The muon beams have the 
particularity that the muon spin and charge are correlated, enabling measurements of the 
DVCST-BH interference via correlated beam charge-spin asymmetries. In addition to the 
COMPASS spectrometer, exclusivity will be determined by detecting the recoil protons in a 
scintillation array surrounding the target. The expected (correlated) range for DVCS and 
exclusive vector meson production is xbj G (0.03, 0.25) and G (1.5, 7.5) GeV^. A future 
electron ion collider, with luminosity several orders of magnitude higher than HERA would 
greatly expand the reach of GPD studies. Maximizing the luminosity is essential to measure fully 
differential cross sections in all kinematic variables. A collider can deliver both longitudinally 
and transversely polarized beams without the accompanying background of unpolarized nuclei 
of polarized targets. A collider would also offer enhanced opportunities for spectator tagging to 
measure neutron GPDs, and recoil tagging for nuclear GPDs. 



6.3. Conclusions 

Deep virtual exclusive scattering offers the tantalizing prospect of forming spatial images of 
quarks and gluons in the nucleon. The GPD formalism has already given us new insight into 
nucleon structure, with evidence for quark angular momentum emerging from GPD models and 
lattice calculations, and global analysis of forward parton distributions and electromagnetic form 
factors. A very important study of DVCS and DES in the valence region has started with JLab 
at 6 GeV and will expand with the 12 GeV upgrade. Several systematic analysis demonstrate 
the important constraints on individual GPDs of the proton and neutron obtained from the 
data [m [THl |T06] . The unprecedented quality of the CEBAF continuous wave beam is essential 
to achieving full exclusivity at high luminosity. The revolution in polarized beams and targets 
over the past two decades allows us a full study of the spin degrees of freedom of DES. Over the 
next decade, JLab and COMPASS will obtain new precision DVCS data spanning a factor of 
20 in XBj, and at each value of rcBj, a factor of two in Q^, with maximal from 4 to 10 GeV^. 
The present JLab data are fully differential in Q^, xbj and t, allowing a systematic study of the 
approach to scaling in both cross section and asymmetry observables. 
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